Quantum Properties of a Nanomechanical Oscillator 
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We study the quantum properties of a nanomechanical oscillator via the squeezing of the oscillator 
amplitude. The static longitudinal compressive force Fo close to a critical value at the Euler buckling 
instability leads to an anharmonic term in the Hamiltonian and thus the squeezing properties of 
the nanomechanical oscillator are to be obtained from the Hamiltonian of the form H = a^a + 
/3(a^ + a)*/4. This Hamiltonian has no exact solution unlike the other known models of nonlinear 
interactions of the forms , {a^a)^ and a^''' + a''' — {a^^ + a^^) previously employed in quantum 
optics to study squeezing. Here we solve the Schrodinger equation numerically and show that in- 
phase quadrature gets squeezed for both ground state and coherent states. The squeezing can be 
controlled by bringing Fq close to or far from the critical value Fc- We further study the effect of 
the transverse driving force on the squeezing in nanomechanical oscillator. 
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42.50.Dv, 05.45.-a, 62.30.+d 



I. INTRODUCTION 

There is currently a wide effort to observe quantum 
behavior in nanoscale devices [E S H, Hi- ^^e limit 
of high resonator frequency with high mechanical quality 
factors and long coherence lifetimes, the nanomechani- 
cal oscillator (NMO) phonons will be analogous to pho- 
tons in an electromagnetic cavity. With current technol- 
ogy it is possible to reach resonator frequency of GHz 
order 5]. At a temperature of around 50mK, one can 
principally prepare the resonator into the ground state. 
These sub-Kelvin temperatures are well within the range 
of todays dilution refrigerators. However cooling the res- 
onator down to these temperatures requires some other 
techniques [1, 0|- 

With the assumption that quantum mechanics should 
apply to these mesoscopic systems, variety of methods 
and techniques have been proposed to observe some 
quantum optical effects including solid-state laser cool- 
ing 't'I, quantum nondemolition measurement 
phonon lasing [10] and squeezed state generation [11, 12]. 
There are also proposals analyzing macroscopic quan- 



tum tunneling [l3|, resonant multi-phonon excitations 
[3] and a variety of methods to entangle mechanical res- 
onators with other quantum systems [isl . [l6l [TtI [T§ ] . 



The next question is what could be the best way to 
study quantum properties of a NMO. In line with the 
work in quantum optics on squeezing, we can consider 
studying the squeezed states of the NMO. One pro- 
posal considers modulating the spring constant to pro- 
duce squeezing [ll| as it is known from the earlier work 
[19I ] that any modulation of the frequency of the oscil- 
lator can result in squeezing. Here we adopt a differ- 
ent model. We consider the situation shown schemati- 
cally in the Fig. [T] We show various forces acting on 
a nanobeam structure which is clamped at both ends 
vibrating in the transverse direction. There is a static 
mechanical force, Fq , acting in the longitudinal direction 
and an ac-driving force to excite vibrations in the trans- 
verse direction. The longitudinal force Fq which is close 
to a critical value at the Euler instability gives rise to 
an additional term in the potential energy which is quar- 
tic in the fundamental mode amplitude x. The effective 
Hamiltonian that describes the system would be in the 
form H = /2m -I- muP'x^ -V /3a:^/4. The derivation of 
this nonlinearity in x is given in the next section. Unlike 
the previous work on squeezing Til in a NMO we would 
consider the effect of the nonlinearity in x. Note that the 
nonlinearity can be switched on and off by controlling _Fo . 
We would thus study the quantized behavior of a NMO 
subject to the force Fq. 



^ AC- driving force, F(t) 

1=^ Static, longitudinal compressive force, 

FIG. 1: (Color online) Schematic diagram of the freely sus- 
pended nanomechanical beam of total length L, width w and 
thickness d. The beam is clamped at both ends. A static, 
mechanical force -Fo compresses the beam in longitudinal di- 
rection controlling the nonlinearity. An ac-driving force can 
be used to excite the beam to transverse vibrations. 



The organization of the paper is as follows. The model 
is described in section [III and the effective Hamiltonian is 
derived briefly by referring to the previous works for the 
doubly clamped elastic rectangular beam. We discuss the 
previous works on the squeezing in nonlinear oscillators 
in section IIIIl Then, we study the quantum dynamics 
and analyze the squeezing properties in section |lVl The 
conclusion and future perspective are given in section [V] 
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II. THE MODEL 

We start with an elastic rectangular beam of length 
L, width w and thickness d as shown in Fig. [TJ The 
beam is freely suspended and clamped at both ends. 
The transverse motion in the direction of d is allowed. 
The dimensions are such that {L w > d) there is no 
appreciable vibrations in other directions. A static me- 
chanical force Fq acts on the beam in the longitudinal 




Here p — m/ L is the mass density, /i = EI is the product 
of the elasticity modulus E and the moment of inertia /. 
In (/)', prime denotes partial derivative with respect to s, 
i.e. d(f>/ds. For small oscillations |0'(s)| <C 1, the La- 
grangian is quadratic and it leads to the linear equation 
of motion 

p4> + p(j>"" + Fo0" = 0. (2) 

The equation of motion can be separated and trans- 
formed into an eigenvalue problem with boundary condi- 
tions applied to the endpoints. One can write the gen- 
eral solution as a superposition 4>{s,t) = J2n't'n{s,t) = 
X]„ -4„(t)5„(s), where gn{s) are the normal modes which 
follow as solution of the characteristic equation and n = 1 
is called the fundamental mode. For the doubly clamped 
nanobeam, we have 0(0) — 4){L) — and 0'(O) = 4)'{L) = 
0. The expressions for (/)„'s are then given by a superpo- 
sition of trigonometric and hyperbolic functions, and the 
eigenfrequencies w„'s are the solutions of transcendental 
equations. For the doubly clamped boundary conditions, 
when Fq is close to the critical value Fc = pijijEf' ^ one 
can get the simplified fundamental frequency iO\(FQ 
Fc) — V^ujQ where luq = {2^"^ /3)^/Ejp{d/ L"^) is the fun- 
damental frequency of the relaxed beam(Fo = 0). The 
parameter e = {Fc — Fq)/Fc is called the distance to the 
critical force and the system reaches to the well known 
Euler instability when e <C 1. The dynamics at low en- 
ergies close to the Euler instability will be dominated 
by the fundamental mode alone. The fundamental mode 
gi{s) can also be expanded close to the Euler instabihty 
in zeroth order in e 

ffi(.)^sin^(^). (3) 

Since the fundamental frequency loi vanishes at the criti- 
cal value Fc, one has to include the contributions beyond 
the quadratic terms 0'^ and 0"^ in the Lagrangian. The 



direction (Fq > for compression). An ac-driving field, 

F{t) = f cos{uJext) can also be added to excite the vi- 
brations. The dynamics of the beam can be completely 
described by the transverse deflection 0(s) parametrized 
by the arclength s G [0, L] in a classical picture. Assum- 
ing single transverse degree of freedom for simplicity the 
nonlinear Lagrangian of the system, for arbitrary strong 
deflections 0(s) is then 0,[2ll, 



(1) 



I 

next higher order terms, — (/i/2)0"20'2 and (Fo/4)0'^, are 
quartic in the Lagrangian. Inserting the normal mode 
expansion 4>[s,i) — ^'^^ Lagrangian 

and assuming that the fundamental mode n — 1 dom- 
inates the dynamics (by neglecting the higher modes 
n = 2,3, . . .) one can quantize the theory by introduc- 
ing the canonically conjugate momentum p = —ihd/dAi 
with the "coordinate" x = Ai- Note that when the driv- 
ing frequency is close to the fundamental frequency of the 
beam, the fundamental mode will dominate also in ab- 
sence of a static longitudinal compression force Fq. How- 
ever, a compression force close to a critical value helps to 
enhance the nonlinear effects which are of the importance 
of this paper. 

By using the above definition of coordinate and the 
conjugate momentum, an effective quantum mechanical 
time-dependent Hamiltonian results describing the dy- 
namics of a single quantum particle 

with the effective mass m* — 3pL/&, the fundamental 
frequency wi = y^(2n'^ /3)^/E/p{d/ L"^) and the nonlin- 

earity parameter /3 = (7r/F)'*F;i(l + 3e) [22]. Now, Eq. 

can be put in a second-quantized form by replacing 
X and p with the creation and annihilation operators 
and a, 

^ = \^ + (5) 

V m*wi V2 

p — \/ m*fkui—={a' — a). (6) 

Upon scaling the Hamiltonian by hcoi we obtain the di- 
mensionless form. 
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H{t) = a^a + 



V2 



fcOs{uJext) 



V2 



(7) 



r 



with the redefined dimensionless parameters, 



/ 



huji ' 



(8) 



where a;o — \/hJmFLO\. One can obtain an expression 
for /3 that depends on the dimensions and the material 
properties of the beam. By substituting the parameters 
in Eq. ([8]) one finds 



/3 = 



1 



1 1 + 3e 



2Lw<P 



,3/2 



(9) 



Note that the equation ^ is valid for e ^ 1 and (3 can 
be controlled by fine tuning the distance parameter e at 
this regime. Table |T] lists the range of the nonlinearity 
parameter (3 as well as the relaxed fundamental frequen- 
cies and the critical compressions for three different sizes 
of Si nanobeams. Note that one should have extremely 
precise control over e to increase the nonlinearity. In light 
of the measurements done in the experiment [23| e was 
found to be of the size ~ 10~^ for a 100 nm length carbon 
nanotube. 



TABLE L Table shows calculated (3 values, the nonlinearity 
parameter, for two different beam size. For each size, the 
relaxed fundamental mode frequency uio and the critical force 
Fc defined in the text, are also shown. The parameter e shows 
the distance to the critical force near the Euler instability. 

Si bar, p = 2330 kg/m^ E = 169 GPa 
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III. PREVIOUS WORKS ON THE SQUEEZING 
IN NONLINEAR OSCILLATORS 



The squeezing produced by the nonlinearity has been 
investigated in the past by using a number of approxima- 
tions however none of these are suitable for the problem 
of the NAMO. Milburn dropped all phase sensitive terms 
from ([7]) and studied [13 the simplified Hamiltonian, 



(10) 



By solving exactly the phase space distribution function 
he showed that squeezing can be obtained for a coherent 
state of amplitude a — 0.5 for very short times. Buzek 
[25| and Tanas [26[ studied Hamiltonian models of the 
form. 



H = hhja^a + i?i/3(ata)^. 



(11) 
(12) 



They solved the Heisenberg equations of motion exactly 
and showed periodic squeezing for coherent states in 
both quadratures. Tanas [l^ showed also that maximum 
squeezing can be obtained in the limit of large mean num- 
ber (lap ^ 1) and small times {t ^ 1). No squeezing 
is allowed for the vacuum (ground) state in the above 
models. 

In this paper, we calculate the squeezing in an an- 
harmonic oscillator for the anharmonicity quartic in x 
which is the amplitude of the fundamental mode of the 
oscillation. One could write a;^ in terms of creation and 
annihilation operators as follows: 



cx (a^ + af = flt^ + a^ + f{a\ a) 



(13) 



where f{a\a) is a polynomial in a and of order four 
which is given by 

/(at, a) = 6 a)^a^ + A{a)^a)a + aa^) + 6(0+2 ^ ^2^ 
-t- 12a+a + 3. (14) 

The Hamiltonian containing the second and third terms 
given in Eq. (|13p give rise to two-photon (or phonons in 
quantum mechanical descriptions of solid systems) tran- 
sitions. It is known that two-photon transitions are nec- 
essary for producing squeezed states and thus the terms 
a+aa^, 0? and their hermitian conjugate should be impor- 
tant. In the literature, multi-photon processes have also 
been analyzed to study normal and higher order squeez- 
ing in the limit of small times [27|, I2M [2^ . In relation 
to this paper, Tombesi and Mecozzi [30] studied the har- 
monic oscillator model which has four-photon transitions 
in the interaction term, 



Hi = /3[a+4 + a" - {a 



t2 2 I 2 t2 

' a + a a' 



)]• 



(15) 



This model was solved exactly. They showed that sig- 
nificant amount of normal and higher order squeezing is 
possible for initial coherent states of amplitude |a| > 1 
with certain phases and for short times. No squeezing is 
allowed for the vacuum and fiuctuations diverge as time 
grows. 
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FIG. 2: Squeezing of the ground state in the anharmonic oscillator model given by Eq. ([7]) changing from weak to strong 
nonlinearity j3; (a) /3 = 0.1, (b) /3 = 0.3, (c) /3 = 0.5, (d) /3 = 1.0. Solid line, Sx{i), normally ordered normalized fluctuations 
in X and dashed line, Spit), normally ordered normalized fluctuations in p as given by Eqs. (|18|) and (|19|l respectively. 
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FIG. 3: The squeezing S'^(t) for the coherent state of amplitudes, (a) |a| = 0.1, (b) |a| = 0.5, (c) |a| = 1.0, (d) \a\ = 2.0. The 
phase of a is equal to 7r/2. The nonlinearity parameter (3 is taken as 0.1. 
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In the description of NMO systems the anharmonic 
Hamiltonian models given in Eqs. ^TU\i . (|lip . and 
(fT4| are no good for vacuum (ground state) squeezing. 
Moreover, the nonhnearity A cannot be controlled exter- 
nally since it is an intrinsic property of the medium. To 
observe the nonclassical (quantum) properties of a meso- 
scopic system in general, the control of the parameter 
that gives rise to the nonclassical behavior would be cru- 
cial for an experimentalist. 

The harmonic oscillator having the nonhnearity of Eq. 
(|13p that we shall work in the next sections, gives im- 
portant squeezing in the in-phase quadrature for both 
ground state and coherent states. Furthermore, it will 
be shown in section III that the physical model of the 
NMO allows one to control the nonhnearity by the appli- 
cation of a static external force. The numerical solution 
of the Hamiltonian shows that the ground state squeezing 
displays periodicity and it stays squeezed for the whole 
cycle of the period. In fact, the ground state squeezing is 
important for the mesoscopic resonators because bringing 
the harmonic oscillator representing the nanomechanical 
system to its vibrational ground state is a necessary pre- 
requisite for quantum state engineering. The effect of 
driving term is also examined. 



IV. QUANTUM DYNAMICS AND SQUEEZING 

We first consider the case in which there is no driving. 
The Hamiltonian 



t 1 o / a + 



(16) 



has no analytical solution. For the numerical calcula- 
tions, we employ the split-operator method [3l| for the 
time propagation of the initial state. In this method, one 
can split the propagator on a time step At as 

U(t + At) = e-'(-f^"+^)'^* (17) 

^ ^-^HoAt/2^-^VAt^-^HoAt/2 ^ 0[{At)% 

where Hq = a+a -1- 1/2 and F = /3(at + a)'*/4. That 
means splitting the exponential of the operators which 
are not commuting is accurate to second order in the 
time step At. Therefore one can make the calculation as 
accurate as possible by taking the time step sufficiently 
small. Now we are ready to analyze the dynamics of 
the oscillations of the fundamental mode. We take the 
ground state of the Hamiltonian Hq as our initial state. 
Then we analyze the dynamics by the application of com- 
pressive force Fq suddenly (in a time interval much faster 
than the oscillator's frequency) in the regime near to the 
Euler instability where the nonhnearity parameter p can 
go from negligibly small values (where the dynamics is 
largely governed by ) to appreciable values. One can 
see how the nonhnearity parameter /3 changes by tun- 
ing the distance parameter e in the table [H We examine 
the dynamics by calculating the time evolution of the 



width of initial wave function in the "coordinate" -space, 
i.e. cc-space (we take the amplitude of the transverse de- 
fiection of the NMO as our coordinate as explained in 
section II.) and the corresponding momentum-space, i.e. 
p-space. Next, we define the squeezing factors 



S^t) - 
Spit) = 



Ax{t) 
Axq 

Apo ' 



(18) 
(19) 



where Ax{t) = {{x'^{t)) - {xit))'^)^/'^ , and Ap{t) = 
((p2(i)) - (p(t))2)i/2, Aa;o and Apo are the widths in 
coordinate and momentum space of the initial wave func- 
tion. The squeezing occurs when one of the expressions 
in Eqs. [18] and [19] get less than one and it is said to 
be perfect when it gets zero. Fig. [2] shows squeez- 
ing for the ground state for different nonhnearity values 
(3 = 0.1,0.3,0.5 and 1.0. 




FIG. 4: Maximum squeezing in x vs. coherent state amplitude 
\a\. The dependence is shown for different values of the phase. 
The nonhnearity parameter /3 is equal to 0.1. 
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FIG. 5: Time evolution of the uncertainty product 
Aa:(t)Ap(f)/Aa;oApo for the ground state wave function of 
7^0 under the dynamics of the Hamiltonian H given in Eq. 
(HI. 



One could also start with a coherent state to check 
squeezing. To prepare the coherent state, just as in 
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FIG. 6: The effect of the driving term on squeezing, (a) shows the squeezing without driving and (b), (c), (d) show the effect 
for the dimensionless driving parameter (/) values of 0.1, 0.2 and 0.3 respectively. The nonlinearity parameter /3 is equal to 
0.1. Solid line is for Sx{t) and the dashed line is for Sp{t). 



the preparation of ground state mentioned in the previ- 
ous paragraph, one can first start with the linear har- 
monic oscillator Hamiltonian Hq. Then the coherent 
state can be obtained by resonantly driving (e.g. by a 
harmonic electromotive force) (3^ the systems ground 
state. Thereafter, by switching the nonlinearity suddenly 
the dynamics of squeezing can be obtained. We ana- 
lyzed squeezing for different initial coherent state ampli- 
tudes and phases, a = laje*^, at the nonlinearity value of 
(3 = 0.1. Fig. [3] shows the time evolution of the squeez- 
ing for increasing coherent state amplitude a at the phase 
— 7r/2. Fig. [4] shows the dependence of the maximum 
squeezing to the phase. It can be seen from the figure 
that the phase does not make much difference for low 
amplitudes (|a| < 0.5) but it changes squeezing behavior 
drastically for amplitudes larger than 0.5. Increasing the 
amplitude increases the maximum squeezing whereas the 
state never becomes squeezed after a short duration. On 
the other hand, low amplitudes shows periodic squeezing 
all the time at a moderate value. 

One can also plot the time evolution of the uncertainty 
product Ax{t)Ap{t) / AxoApo for the ground state wave 
fimction of the linear Hamiltonian Hq. Fig. [5] shows that 
the oscillator recovers its minimum uncertainty periodi- 
cally and the fluctuation remains bounded. 

Next we analyze the effect of driving in the dynam- 
ics of the Hamiltonian given in Eq. (O. We calculate 



the propagator again by using the split-operator method. 
This time, we include the nonlinear term into Hq and we 
take the time dependent driving term as V{t) to employ 
the splitting given in Eq. p?)) . Fig. [6] shows the time 
evolution of the squeezing factors given in Eqs. (fT5|) and 
for the in-phase and the out-of-phase quadratures 
for driving parameter values of / = 0.0, 0.1, 0.2 and 0.3 
at the nonlinearity /? = 0.1. We take the ground state 
of the linear Hamiltonian as the initial state. The fre- 
quency of the driving term, LUex , is on resonance with the 
frequency of the oscillator, uji. As clearly seen, the effect 
of the driving is to enhance the squeezing in x periodically 
to a larger value. For the values of / in consideration, 
maximum squeezing goes from 17% (for / = 0.0) to 30% 
(for / = 0.3). Increasing the driving further does not 
enhance squeezing, rather it starts to enhance the fluc- 
tuations and the squeezing cycle becomes shorter. And 
the driving frequencies off the resonance does not help 
enhancing the squeezing at all. The upper limit fmax 
for the regime of validity of the fundamental mode de- 
scription is given by the first harmonic threshold, that is 
fmax < ^2 ~ Shuji [13] which means for the dimension- 
less parameter, fmax < 3. 

We include the damping effects in the master equation 
formalism of the form 

P = -i[H, p] + -f{apa'^ - ^ a^'ap - ^ pa^a), (20) 
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FIG. 7: (Color online) Squeezing in the driven damped anharmonic nanomechanical oscillator. The red and blue lines show 
S-p{t) and Sxit) respectively, (a) and (b) show squeezing for the ground state with the damping constants 7 = O.OOIli^i and 
7 = O.OOSoJi respectively, (c) and (d) show squeezing for the thermal state with the damping constant 7 = O.OOScJi at the 
temperatures kT = O.lhoJo and kT = 0.5huJo respectively. The other parameters are kept at /3 = 0.1 and / = 0.1 for all. 



where H is the Hamiltonian given by the Eq. ([7]) and 
7 is the damping constant. Here, 7 is the rate for the 
oscillator to decay from the first excited state to the 
ground state. The rate for the decay from level n to 
n — 1 is n^. After solving the master equation numer- 
ically we plot the squeezing factors Sx{t) and Sp{t) in 
Figs. [7]Ja) and[7]Jb) for the ground state. As seen clearly, 
the damping effect shows itself as reducing the fluctua- 
tions exponentially to a steady state value. The master 
equation formalism also allows us to analyze the squeez- 
ing for the thermal state of the relaxed beam (described 



by the linear Hamiltonian Hq), p{t — 0) 



-Ho/d 



where Ho = afa + 1/2, 9 = kT/fuvo and Z = rr(e~^°/^). 
Figs. [7[c) and Uld) show that one need not to start at 
K to obtain squeezing. As long as one keeps the noise 
of the environment sufficiently low, at higher tempera- 
tures of the beam {kT of fractions of the level spacing 
fkoo), although there is no substantial squeezing initially, 
at steady state, the amount of "coordinate" squeezing is 
reaching to the ground state value. 



V. CONCLUSION 

In this paper we show that squeezed states can be ob- 
tained in the amplitude of the fundamental mode of a 



nanomechanical oscillator which has quartic nonlinear- 
ity in its effective single particle quantum mechanical 
Hamiltonian. The quantum dynamics of the system is 
solved numerically both with and without external ac- 
driving. For various strength of nonlinearity and driv- 
ing, the squeezing dynamics is investigated for both ini- 
tial ground state and coherent states of different ampli- 
tudes. The terms that lead to multilevel transitions in 
the quartic nonlinearity is compared to similar models 
and proved advantageous results. It should be noted 
that working closer to Euler buckling instability produces 
larger squeezing due to large nonlinearity. This is remi- 
niscent of similar result for quantum optical systems [ssj . 

Finally, we comment on the possibility of observing 
squeezing in nanomechanical beams. The obvious way to 
detect squeezing seems to be one of the displacement de- 
tection methods. There are different displacement sens- 
ing techniques such as optical, magnetomotive and single- 
electron transistor (SET) technique. Among these read- 
out strategies, the SET -demonstrated by many research 
groups to be a very sensitive detector of charge- enables 
position detection for NMO devices in their low energy 
quantum states. In this detection method, a charged 
NMO is capacitively coupled to an SET. The oscilla- 
tor's motion induces a change in the charge on the gate 
electrode of the SET and the change in the SET's con- 
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ductance can be directly monitored. Demonstrations of 
position detection of NMO's are developing. The first 
experiment [s^ reached the position sensitivity within a 
factor of 100 of the quantum limit IS.xsql — \/h/2mujQ. 
A year after another group improved this sensitivity 
to a factor of 4.3 of the quantum limit set by the Heisen- 
berg uncertainty principle of a Nanomechanical resonator 
nearly its ground state. These demonstrations suggests 
that beating the standard quantum limit in the position 
detection technology is highly probable in future. Once 
this is achieved, by reading out the position each time 
with a different nonlinearity parameter one can prove the 
existence of tunable squeezing. There are also proposals 
suggesting the positioning of a NMO inside an electro- 
magnetic cavity. It is proposed that the interaction of 
cavity photons with a NMO in a squeezed state would 



produce a nonclassical statistics of photons at the exit. 
So one can easily detect squeezing simply by quantum 
statistical analysis of photons at the output of the cav- 
ity. 

In conclusion, in nanomechanical beams described in 
this paper, one can remarkably control the quantum 
squeezing externally just by controlling the strain pro- 
vided by a classical compressive force. However, to con- 
trol the nonlinearity at the desired values, one would 
have to apply compression with extreme delicacy, Fc — 
Fq ~ 10~^Fc. Controlling the strain to this precision 
for sufficient time to identify squeezing may be diffi- 
cult. Thus, while observing squeezing will be challenging, 
the prospect of exploring tunable quantum squeezing in 
nanomechanical beams and the connection to Euler buck- 
ling instability are intriguing. 
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